Aperiodic Ising model on the Bethe lattice: Exact results 
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We consider the Ising model on the Bethe lattice with aperiodic modulation of the couplings, which 
has been studied numerically in Phys. Rev. E 77, 041113 (2008). Here we present a relevance- 
irrelevance criterion and solve the critical behavior exactly for marginal aperiodic sequences. We 
present analytical formulas for the continuously varying critical exponents and discuss a relationship 
with the (surface) critical behavior of the aperiodic quantum Ising chain. 

PACS numbers: 64.60.-i, 64.10.+h, 05.50.+q 



I. INTRODUCTION 

Disorder and different types of inhomogeneities are in- 
evitable features of real materials. Their presence may 
modify the physical properties of a system and their effect 
can be particularly strongclose to singularities, such as 
at phase transition points [HQ. -"- n some cases the pertur- 
bation can change the universality class of a second-order 
phase transition. In this respect relevance or irrelevance 
of an inhomogeneous perturbation can be analyzed in 
terms of linear stability at the pure system's fixed point 
as first performed by Harris [3( for uncorrelated bond 
disorder. The classification of the critical behaviors of 
disordered systems with a random fixed point is a chal- 
lenging and theoretically very difficult task. 

Another type of inhomogeneities is introduced by 
quasiperiodic or, more generally, aperiodic modulations 
of the couplings. Since the discovery of quasicrystals [4] 
and due to the progress in molecular beam epitaxy, allow- 
ing for the preparation of good quality multilayers with 
a prescribed aperiodic structure [5|, there has been an 
increased interest to study theoretically the phase tran- 
sitions in such nonperiodic systems (fj. These systems 
can be considered as somehow intermediate between pure 
and random ones and are expected to display a rich va- 
riety of critical behaviors. Indeed a generalization of the 
Harris criterion predicts that, depending on the strength 
of the fluctuations of the aperiodic sequence and the 
value of the correlation length critical exponent v of the 
pure system, an aperiodic perturbation may be irrele- 
vant, marginal, or relevant 0, A series of works on 
the critical behavior in different aperiodic systems con- 
firms the validity of the _generalized relevance-irrelevance 

criterion 0, M M M M El, EE El, HE El EI H IH • 

Interestingly in the presence of aperiodicity one can 
observe truly marginal behavior. Then the critical ex- 
ponents are nonuniversal and their value varies contin- 
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uously with the amplitude of the aperiodic perturba- 
tion. Such a behavior has been obtained exactly for 
the aperiodic Ising quantum chains using a renormaliza- 
tion group transformation or a finite-size scaling analysis 
0, El Ell, El, El, HI, as well as for the interface delo- 
calization transition in the Penrose quasiperiodic lattice 
[HI]. Nonuniversal critical behavior is expected to oc- 
cur in real higher-dimensional systems, too, for example, 
in a three-dimensional tricritical system where v = 1/2. 
In this case, however, no exact results are available yet. 
Numerical studies of the related mean-field model with 
Fibonacci modulation of the couplings show nonuniversal 
critical behavior [H [HI . 

More recently the Bethe-lattice Ising model, which also 
belongs to the mean-field universality class [24j , has been 
studied numerically for two types of perturbations [25| . 
With a Fibonacci modulation of the couplings, the clas- 
sical mean-field exponents are recovered whereas, for a 
period-doubling (PD) modulation, the magnetic expo- 
nents are nonuniversal. The difference in the relevance 
of the Fibonacci modulation for the two mean-field mod- 
els is due to the different ways in which the mean-field 
behavior is realized, but this question has not been stud- 
ied so far. 

In this paper we continue the study of the aperiodic 
Bethe-lattice Ising model. Our motivations are twofold: 
First, we are interested in the formulation of a relevance- 
irrelevance criterion adapted to this system in order 
to explain the conflicting results of previous numerical 
works. Second, and our more important motivation, we 
can provide an exact solution of the problem and in this 
way we obtain analytical formulas for the continuously 
varying critical exponents, among others for the PD se- 
quence studied before numerically. 

The structure of the paper is the following. The model, 
the aperiodic sequences and the corresponding relevance- 
irrelevance criterion are presented in Sec. [TTJ The crit- 
ical behavior of marginal aperiodic sequences is studied 
analytically in Sec. IIII1 in which we point out a close 
relationship with the (surface) critical properties of ape- 
riodic quantum Ising chains. Our results are discussed 
in Sec. IIVI and details about the calculation of sums of 
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unusual [26|: There is no long-range order at T > and 
the free energy displays a power-law singularity in H with 
a critical exponent increasing smoothly from 1 to infinity 
as T goes from to the Bethe-Peierls temperature. 

The Bethe lattice may be embedded in a d = oo hyper- 
cubic lattice by placing each link into a different lattice 
direction. Then the geometrical distance L, measured 
in lattice parameter units, between two spins located N 
layers apart grows as 
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N 



(4) 



FIG. 1: Aperiodic Bethe lattice with coordination number 
z = q + 1 = 3. The couplings K n /q are aperiodically modu- 
lated. 



since the L steps are mutually orthogonal. 



B. Aperiodic sequences 



aperiodic variables are given in the Appendix. 



II. APERIODIC PERTURBATION AND ITS 
RELEVANCE 

A. Hamiltonian 

We consider the spin 1/2 Ising model on a Bethe lattice 
with coordination number z = q + 1. The Hamiltonian 
reads 
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n (j,i)ei 



Kn 
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(i) 



where f3 = is the inverse temperature. The first 

sum runs over the successive layers of the lattice indexed 
by n (see Fig. [TJ and the second over the bonds between 
the sites in layers n + 1 and n. The couplings K n /q 
are aperiodically modulated and properly normalized in 
order to allow us to take the mean- field limit q — > oo. 
They are parametrized as 



Kn 



fn = 0, 1 



(2) 



where r is the ratio of perturbed to unperturbed cou- 
plings and the binary variables /„ follow some aperiodic 
sequence. 

Let O n = (Oi) be the thermal average of some local 
operator Oi in layer n. This average is fluctuating from 
layer to layer due to the aperiodic modulation of the cou- 
plings. Following Ref. [25|], a mean value O is defined by 
giving the same weight to the different layers, 



1 N 



(3) 



n=l 



With this choice, the usual mean-field critical behavior is 
obtained for the unperturbed system. When each layer 
is weighted according to its number of spins, the surface 
spins dominate the critical behavior which is then quite 



As in Ref. [25| we consider aperiodic sequences gener- 
ated via substitutions on the binary digits and 1. For 
the PD sequence [1| we have 5(0) = 1 and 5(1) =00 
which, starting on 0, give successively 



0. 

1, 

10 0, 

1 1 1. 



(5) 



The properties of the sequence can be deduced from its 
substitution matrix [29|, [3(| with entries Mij giving the 

numbers nf of digits of type i in S(j). In the case of 
the PD sequence one obtains 



M 



5(0) 5(1) 
5(0) 5(1) 



11 



1 2 
1 



(6) 



The entries in M p give the numbers of digits of each type 
in the sequence after p iterations. The length N of the 
sequence obtained after p iterations (which is also the 
number of layers on the Bethe lattice) is related to the 
leading eigenvalue Ai of the substitution matrix through 
N = A[. Let 



n N 



N 

E 

71=1 



fn , Poo = Jam , (7) 

N—>oo IV 



be the number of 1 in a sequence with iV digits and the 
corresponding asymptotic density, respectively. On the 
Bethe lattice, according to Eq. ([3]), the mean value of 
the coupling is given by 



— 1 N 

K = lim — V K n = K+ lim — K(r 



N^OG N 

n=l 
K + Poo A, 
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where A = K(r — 1) is the amplitude of the aperiodic 
modulation of the couplings. The mean deviation from 
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if on a system with N ~ layers takes the form 
1 N 



N 

— K% ~ A TV" 

TV 2 



(9) 



where u) is the wandering exponent of the sequence given 
by 



taJAal 
InAi 



(10) 



in terms of the second leading eigenvalue A2 of the sub- 
stitution matrix. For the PD sequence, according to Eq. 
(O, we have Ai = 2, A 2 = — 1, so that uj = 0. 



C. Relevance-irrelevance criterion 

The Harris argument [3j, showing that thermal ran- 
domness is a relevant perturbation only when the spe- 
cific heat exponent a of the pure system is positive, has 
been generalized to the case of aperiodic perturbations 
in Refs. 0,0. 

The argument can be adapted to our problem as fol- 
lows: Near the critical point of the pure system, the rel- 
evant length is the correlation length £ diverging as t~ v 
where t ~ \K — K c \ measures the deviation from the 
critical temperature. To the length £ is associated the 
number of layers N = £ 2 according to Eq. Q and the 
aperiodic perturbation induces a shift in the critical tem- 
perature St ~ SK(N) which, according to Eq. ([9]), takes 
the form Si ~ f 2 ^- 1 ) ~ The ratio 



A continuous variation of the magnetic exponents was 
observed in Ref. [25[ for the PD sequence with ui = 0. 
This marginal behavior is expected since yt — 1 jv = 2 
for the mean-field Ising model. On the contrary, the Fi- 
bonacci sequence, with u> = — 1 leads to an irrelevant 
perturbation. It does not change the critical behavior 
which remains classical. 



III. CRITICAL BEHAVIOR 

A. Finite-size behavior of the magnetization 

We consider an n + 1-generation branch defined as an 
initial site with spin cr n +i connected to q = z — 1 re- 
generation branches as shown in Fig. [TJ a one-generation 
branch is a single site. Let Z„ be the sum of the con- 
tributions to the partition function of an n-generation 
branch with initial spin either up (+) or down (— ). It 
satisfies the recursion relation 



Z± +1 = e ±H (e±«^Z+ + e^Z n 



(14) 



In the mean- field limit q — ► 00, the nth layer magnetiza- 
tion m n — (<r n ) may be written as 



Zn + Z n 



(15) 



since the contribution of the single branch going forward 
can be neglected compared to the contributions of the q 
branches going backward. Expanding the exponentials 
in Eq. (fl~4"l) one obtains 



St 

T 
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1 + 2z/(w - 1) 



(11) 



gives the relative strength of the aperiodic perturbation. 
It diverges, and thus the perturbation is relevant, when 
the crossover exponent <f> > 0. It is irrelevant when <fi < 
and marginal when (f> — 0. In this latter case the ape- 
riodicity may lead to a nonuniversal behavior with some 
exponents varying continuously with the amplitude of the 
perturbation. 

The same result can be obtain by studying the scaling 
behavior of the perturbation amplitude A in Eq. 
Under a change of the length scale by a factor b = L/L', 
Eq. (gl) leads to N' = N/b 2 and 5K(N), with scaling 
dimension y t = and transforms as 



(SK)' ~ A' Af" 



bfAN" 



so that 



A' = 6f f + 2 (" _1 ) A . 



(12) 



(13) 



Thus the scaling dimension of A is 4>/v and the per- 
turbation grows under rescaling (is relevant) when <f> is 
positive. 
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(17) 



and the layer magnetization satisfies the recursion rela- 
tion 



m n+ i = tanh(ii + K n m n ) . 



(18) 



Expanding to the first order in the external field H and 
to the third order in the first layer magnetization mi, 
one has 



m n+ i = HjJjf^JJirri+miJJjri 

i=l ft=lj=l i=l 
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i=i fe=ij=i 



(19) 
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According to Eq. ([3]) the mean value of the magnetization 
on a system with n layers is given by 
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(20) 
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(21) 



The critical point of the system is obtained by analyz- 
ing the asymptotic behavior of 01. It is divergent (goes 

r n \ 1/N 

to zero), if (IX^i^ ) i s greater (smaller) than 1. 
Consequently the critical point is given by the condition 



1 N 

lim — V ]nKi ■ = In if = , 



leading to 
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(23) 



B. Relation with the one-dimensional Ising model 
in a transverse field 

Let us consider the inhomogeneous quantum Ising 
chain with Hamiltonian 



1 



1 1 



(24) 



where sf and sf are the components of a Pauli spin opera- 
tor associated with site I, Ji is the first-neighbor exchange 
interaction, and h is the transverse field. On a chain with 
size L and the end spin fixed, the surface magnetization 
m s satisfies the relation 1311 



= 1 



En* 

1=1 i=i 



(25) 



where Ai = Ji/h. This is just the form of the sum giving 
TVai in Eq. (l2Tj) . At the critical coupling, generally given 
by In A = [32], mj 2 scales like L 2Xs where x s — fts/v is 
the scaling dimension of the surface magnetization with 
x s = 1/2 for the unperturbed quantum Ising chain. 

This exponent has been determined analytically in the 
case of marginal aperiodic modulations of the couplings 
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FIG. 2: (Color online) Variation of the critical exponent /3 
of the spontaneous magnetization with the coupling ratio r 
for the PD, paper-folding (PF), and three-folding (TF) se- 
quences. The exponent is minimum and takes its mean-field 
value (3 = 1/2 for the unperturbed system at r = 1. 



EE E3, GH. Since N is replaced by L, 2v = 1 is re- 
placed by fi s i n g = 1 in the expression ([Til of the crossover 
exponent <j). Thus marginal behavior is obtained for the 
same value u> — of the wandering exponent. 
With the parametrization A^ = Ar^ , one obtains 



x s (r) 



hx(r 



1/3 



-1/31 



21n2 



(26) 



for the PD sequence. Changing r into r _1 does not 
change x s for this sequence, but this is not generally true. 

The leading behaviors of the sums appearing in Eq. 
(PH) are calculated in the Appendix. One may notice 
that since Ki corresponds to A~ 2 , r in Eq. ([26]) has to 
be replaced by r~ x / 2 in the scaling exponent of 01. At 
the critical point, the different coefficients scale with N 
as follows 



= 2z s (r 1 / 2 ) + 2z s (r- 1 / 2 )-l, 
= 2x s (r-^ 2 )-l, 

= 2x s (r- 1 ) + 2x s {r- 1 ' 2 )-l. (27) 



C. Finite-size scaling and critical exponents 

In a finite-size system with N layers the free energy 
density g is a function of the deviation from the critical 
temperature t and the external field H, of the system size 
N and also of the magnetization mi of the first layer. 
Under a change of the length scale by a factor b, these 
variables transforms as 



do 


- N X0 , 


Xo 


a\ 


~ N Xl , 




"3 


~ N X3 , 


X3 



t' = b yt t, H' = b VH H, N' 



N 



b 2T m l . (28) 



For a truly marginal system, the thermal dimension y t = 
\/v has to keep its unperturbed value y t = 2. 

The critical exponents f3, governing the temperature 
dependence of the spontaneous magnetization, and 7, 
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FIG. 3: (Color online) Variation of the critical exponent 7 of 
the susceptibility with the coupling ratio r for the PD, PF, 
and TF sequences. The exponent is minimum and takes its 
mean-field value 7 = 1 for the unperturbed system at r = 1. 



governing the temperature dependence of the suscepti- 
bility, have the following expressions 

a d c -y H 2y H - d c 
= — ~ , 7 = . Vt = 2 , (29) 



Vi 



Vt 



where d c is the upper critical dimension of the problem. 
The dimensions yu, t, and d c remain to be determined 
for the aperiodic system. 

The free energy density transforms as 

g{t, H, AT, mi) = b'^g \b v H, b y »H, ^b 2T m^j . (30) 

The mean value of the magnetization m — dg/dH reads 

m(t, H, N, mi) = b- d - +yH m (b Vt t, b VH H, ^,b 2r m^j 



It can be expanded in powers of mi as 
m(t,H,N, mi ) = ^2b yH - d " +2kT 



(31) 



k>Q 



<m (fc) (b yt t,b yH H,?pJm'l.{32) 



At the critical temperature, the leading dependence on a 
small external field H comes from the term of order zero 
in the expansion (|32p . With t = and b 2 — N we have 



m{0,H,N, mi ) = N { y H - d " )/2 ip% ] (N Vh/2 H) + O(mi) 
~ N^y"- d ^' 2 H. (33) 

When H = 0, Yn, is odd in mi so that, with b 2 — N, one 
obtains 

m(t, 0, N, mi) = iV T +(^- rf c)/2^(i) 

+A r3^+to«-^)/2 ¥ ,(3)( iv ^ TO 3 + ... ( 34 ) 



TABLE I: Critical behavior of the Bethe lattice with an ape- 
riodic modulation following the PD sequence: comparison of 
the exact values of the critical exponents j3, 7 in Eq. (|35[) 
and S deduced from the Widom scaling law to the numerical 
values obtained in Ref. 12511. 



r = 1 



r = 2 



r = 7 



/3 (exact) 

/3 (numerical) 

7 (exact) 

7 (numerical) 

S (exact) 

S (numerical) 



1/2 
1 

3 



0.5094795 

0.5093(4) 

1.0192114 

1.0197(2) 

3.0004955 

3.0006(2) 



0.5675859 
0.5664(5) 
1.1491583 
1.1499(4) 
3.0246421 
3.0266(9) 



Comparing Eqs. (J33J) and (34]) to Eqs. (20]) and (27j 
one can deduce the values of the critical exponents 

7 = x Q = 2x s (r 1 / 2 ) + 2x s (r~ 1 / 2 ) - 1, 



/3 = 



2 

x 3 - .Xl 



x s (r- 1 )-2 : r s (r- 1 / 2 ) + l, 
Mr- 1 ), 



2 + 4x s (r" 1 ) + 4x s (r 1/2 ) - 4a; s (r" 1/2 ) . (35) 



For the unperturbed system (r = 1, x 8 = 1/2) the mean- 
field Ising values, 7 = 1, /3 = 1/2, d c = 4, are recovered. 
The variations of j3 and 7 with r for the PD sequence are 
shown in Figs. [2] and [3] Similar results for the paper- 
folding (PF) and three-folding (TF) seq uences are also 
shown [33[. In these cases the values [l7l [l8j 



x a (r) 



ln(l 



2 In 2 



(PF), x s {r) 



In(2 + r) 
21n3 



(TF). 
(36) 

have been used in Eq. ([33]) . 

For the PD sequence, the values of f3, 7, and 6 (deduced 
from the Widom scaling law 5 = 1 + 7 IP) given in Table 
J] for r = 2 and r = 7 compare well with the numerical 
values obtained in Ref. [25| (notice that r in the present 
work corresponds to r — 1 in (25[). 



IV. DISCUSSION 

In this paper we have studied the critical behavior of 
the Ising model on a Bethe lattice with aperiodic modu- 
lation of the couplings. Our first result is the relevance- 
irrelevance criterion of Eq. (TTJ) which is adapted to the 
Bethe-lattice problem. As for the aperiodic quantum 
Ising chain or previous mean- field models, it has a form 
typical of a one-dimensional aperiodicity. The difference 
lies in the fact that the length of the sequence does not 
scale here like a physical length L but like the number 
of layers N on the Bethe lattice, which itself scales like 
L 2 . As a consequence, the correlation length exponent 
v is replaced by 2za Thus for the Bethe-lattice problem 
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with v = 1/2 the aperiodicity is marginal when the wan- 
dering exponent u) = as for the Ising quantum chain 
with ^i s ing = 1. For the same reason, the aperiodicity 
is irrelevant on the Bethe lattice for the Fibonacci se- 
quence with cu = — 1 [25l ] whereas it is marginal for the 
one-dimensional mean- field models [22L [23| . 

We have solved the critical properties of the aperiodic 
Bethe-lattice Ising model exactly and we have observed 
further similarities with the quantum Ising chain. For 
marginal aperiodic sequences, such as the PD sequence, 
the critical exponents are nonuniversal in both cases and 
the Bethe-lattice exponents can be expressed in terms of 
the surface magnetization exponent of the quantum Ising 
chain, taken at different values of the aperiodic coupling 
ratio r. Since for the quantum Ising chain there is a vast 
literature about exact solutions for different marginal se- 
quences, from these we can immediately translate the 
corresponding analytical results for the Bethe lattice. 

We have also noticed that for the Bethe lattice, in order 
to satisfy the scaling relations, a varying upper critical 
dimension d c (r) has to be introduced. As a matter of 
fact, this result follows from an analysis of the Ginzburg 
criterion [34j . too. Something similar occurs for the Ising 
quantum chain, in which th e dy namical exponent z was 
found to be r dependent [l||[l74- As a consequence, here 
also the effective dimension of the system d — 1 + z varies 
continuously with r. 

Our investigation could be extended into several di- 
rections. For relevant aperiodic sequences, such as 
the Rudin- Shapiro sequence, first-order transition is ex- 
pected in one range of the ratio, r < 1, whereas in the 
other range, r > 1, the magnetization should exhibit an 
essential singularity at the critical point [351 ] . instead of 
a power law observed for marginal perturbations. This 
type of essential singular behavior is probably the rule 
for random interactions, too. 
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APPENDIX: CALCULATION OF THE SUMS 
THROUGH RENORMALIZATION 

Let us consider the sum 

N p N 

S N (K, r) = 1 + Krh = J2 KPrHp ' n ° = ' 

p— 1 i—1 p—0 

(A.l) 

such that a\ in Eq. (f2"Tj) is given by Sff-i(K,r)/N. For 
the PD sequence the following relations are satisfied (see 



Ref. [14| where and 1 are exchanged) 
fzk = 1 - fk ) fak+i = , n 2 k = n 2 k+i =k-n k . (A. 2) 

Splitting the sum into even and odd parts and ignoring 
minor end corrections, one obtains 



N/2 



N/2 



S N (K,r) = Y, K 



2k n 2 k 



K 2k+l r n 2k + 1 



k=0 



k=0 

N/2 N/2 

= Y^{K 2 r) k r- nk + K^2(K 2 r) k r- nk 

k=0 fc=0 

= (l + ^S^^r.r- 1 ), (A.3) 

A second iteration leads to 

S N (K,r) = (1 + K){l + K 2 r)S N/i {K A r,r) (A.4) 

which is a renormalization transformation leaving r in- 
variant, dividing TV by 4 and changing K into K' — K 4 r. 
In the infinite system this transformation has a nontrivial 
fixed point at 



K* = K r = r~ 1/3 



(A.5) 



which is the critical coupling of the problem, in agree- 
ment with the general result of Eq. (|23[) . At the critical 
point, Eq. (A4) gives the finite-size scaling relation 

S N (K c ,r) = (l + r-^Xl + r^S^K^r) 

= (r'/' + r-^JV^-r). (A6) 

Injecting the power law Sn{K c ,t) ~ AN 10 ^ into Eq. 
(A6), one finally obtains 



in agreement with the value of x\ given in Eq. (|2T|) . 

The values of Nao and Na,3 in Eq. (f2T|) are given by 
the sum 



T N {K, r;m)=J2 K "" r 
p=i 1=1 



ml mni 



(A. 



with to = — 1 for ao and m — 2 for 03. The leading 
contribution to this sum for large N values is given by 



JV 



T N (K,r;m) ~ ^ K p r n "S p (K m , r m ) 

p=0 



N 



p=0 



(A.9) 



Proceeding as before for Sn, the sum can be split into 
even and odd parts and after two iterations one obtains 

M 

T N (K C , r; to) ~ ^^{r 1 ^ + r^ e )T N/i (K c , r; to) . 

(A.10) 
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It follows that, at the critical point, this sum scales with 
N as 

T N (K c ,r;m) ~ N u{r ' m \ w(r;m) = w(r) + w(r m ) 
w(r;m) = 2a; s (r- 1 / 2 ) + 2x ;5 (r- m / 2 ) . (A.ll) 



With the appropriate values of m one recovers the expo- 
nents in Eq. (|27|) . 
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